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One of the most important developments in condensed matter physics 
in recent years has been the discovery and characterization of graphene. A 
two-dimensional layer of Carbon arranged in a hexagonal lattice, graphene 
exhibits many interesting electronic properties, most notably that the low 
energy excitations behave as massless Dirac fermions. These excitations 
interact strongly via the Coulomb interaction and thus non-perturbative 
methods are necessary. Using methods borrowed from lattice QCD, we 
study the graphene effective theory in the presence of an external mag¬ 
netic held. Graphene, along with other (2 -|- l)-dimensional held theories, 
has been predicted to undergo spontaneous breaking of havor symmetry 
including the formation of a gap as a result of the external magnetic held. 
This phenomenon is known as magnetic catalysis. Our study investigates 
magnetic catalysis using a fully non-perturbative approach. 
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1 Introduction 

The field of graphene has seen tremendous growth since its experimental verification 
in 2006 [Tj. The unusual electronic properties exhibited by monolayer graphene and 
its related structures have been partly responsible for driving this growth [2]. 

What makes graphene’s electronic properties so fascinating is that the electrons 
interact strongly. The graphene fine-structure constant, ag, analogous to the fine- 
structure constant of QED, is of the order of one. Although perturbative approaches 
have had success in elucidating aspects of graphene, one would also like a non- 
perturbative approach in order to investigate previous claims and explore new avenues 
in a controlled, systematic way. We do so here. 

Magnetic catalysis, spontaneous symmetry breaking in the presence of a magnetic 
field, is thought to be a universal mechanism which is particularly applicable to 
field theories in (2 -|- 1) dimensions [3]. Using the low-energy effective field theory 
(EFT) describing graphene, studies applying the Schwinger-Dyson approach in the 
Hartree-Fock approximation find that a dynamical mass is generated. The form of 
the dynamical mass in (2 -|- 1) dimensions is given by rridyn oc ag^^/\eB\ |1]. 

The low-energy EFT has an internal symmetry described by U (4) which is broken 
down to U{2) ® U{2) with the appearance of a nonzero value for the condensate 
{'ijj'ip) and the generation of the dynamical mass. We use a discretized version of the 
EFT which employs staggered lattice fermions to correctly describe Nf = 2 species 
of degenerate four-component Dirac spinors. Initial lattice studies did not see a 
nonvanishing condensate at arbitrary coupling [5]. Through careful attention to the 
zero-temperature and chiral limits, we show in a nonperturbative context that the 
graphene EFT does exhibit magnetic catalysis. 

2 Graphene EFT 

Graphene is predicted to be a semimetal from its band structure. One sees that 
around two inequivalent corners of the first Brillouin zone, called Dirac points, the 
dispersion is linear and the excitations are described by [2] 

'H(fc) = hvp^ ■ k. (1) 

Counting the spin of the electron and the sublattice degree of freedom, one has a 
total of 8 fermionic degrees of freedom which compose the Nf = 2 species of four- 
component Dirac fermions in the EFT. The quasiparticles interact via a Coulomb-type 
interaction described by 

'Hint = j dfj dr'p (r) (A'), (2) 

where p(r) = Tl(r)T(r), Tj = and e refers to the dielec¬ 

tric constant of the substrate on which the graphene sheet is placed. In the definition 


1 




of K and k! refer to the two Dirac points, A and B refer to the two ineqniva- 
lent snblattices, and /x refers to the electron’s spin. The effective conpling between 
qnasip art ides is described by the graphene fine-strncture constant which is given by 
ttg = e^/evpAn > 1. 

3 Results 

We employ one flavor of staggered fermions in (2 + 1) dimensions. Interactions be¬ 
tween the fermions are mediated by a scalar gange potential which lives in (3 -|- 1) 
dimensions. On the lattice this is described by the non-compact U{1) Wilson action. 
For completeness we write onr action below 
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(3) 


(4) 


where i labels the spatial coordinate and we have taken the spatial and temporal 
lattice spacings, and at, to be distinct. The mass term explicitly breaks the sym¬ 
metry and is needed as an infrared regulator. We take the zero mass (chiral) limit to 
recover the EFT. The fermion action in ([^ is known to have 0{a?) splitting in the 
mass spectrum. To suppress this splitting at this order and obtain a better approx¬ 
imation to the continuum action, we replace the links Uo{n) with tadpole-improved 
“fat” links and also add a third-nearest neighbor term. This results in what is known 
as the improved AsqTad action |B]. 

It is known that the EFT undergoes spontaneous symmetry breaking at large 
coupling without an external magnetic field [7] . It was thus necessary Erst to identify 
the symmetric phase as a function of the inverse coupling, (3 = Once the 

symmetric region was identihed, we fixed the coupling within that region. We then 
introduced the external magnetic held at several values of the quantized hux. The 
effect of the external held on the condensate can be seen in LHS Fig. ([^. The 
extrapolation of a to zero as the bare mass is taken to zero signals that the symmetry 
is restored when the explicit breaking term is removed. Typically, in lattice QCD, 
this is due to ehects arising from the hnite spatial size of the box. However, thermal 
ehects can also play a role. The temperature, T = is set by the temporal 

extent of the lattice. As one can see in RHS Fig. ([^, the ehects of a hnite spatial 
volume are not signihcant and one still has a vanishing in the chiral limit. Further 
proof of this came from checking the screening masses of the lightest excitations of 
the theory. 
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Figure 1: (LHS) a = {'ip'ijj) as a function of the bare fermion mass. We show the 
behavior at two values of the magnetic flux, = 0.125 (black points) and = 
0.1875 (blue points) in units of a^, as well as at zero external magnetic held (red 
points). One notices that a vanishes with m in the absence, as well as in the presence 
of, the external held. We argue that vanishing in the presence of an external held is 
a thermal ehect. (RHS) a as a function of mass with varying spatial volume for a 
magnetic hux = 0.125. The lattice volumes are reported in the form xNzX Nt. 

From an analysis of a{T/m), illustrated in LHS Fig. ([^, one can see that hnite 
temperature ehects strongly ahect the chiral extrapolation. Namely, one sees that a 
zero-temperature extrapolation must be performed before the chiral limit is taken as 
hnite temperature ehects become signihcant when T/m > 1. The zero-temperature 
extrapolated points are displayed in RHS Fig. (|^, which lends support that a nonzero 
condensate exists in the chiral limit. 

4 Conclusion and Outlook 

The results of our fully non-perturbative calculation are consistent with magnetic 
catalysis. We have shown that hnite volume and temperature ehects are under control 
and that there is a systematic way of performing both the chiral and zero-temperature 
limits. A nonzero value of the condensate was found at a single magnetic hux after 
performing both of the aforementioned limits. 

Next, we plan to investigate the dependence of the condensate on the external 
magnetic hux and the coupling. This will allow us to compare with the predictions 
of 0. 
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Figure 2: (LHS) a plotted as a function of the scaling variable T/m. At small values 
of T/m, (j tends toward a nonzero value while at large T/m the condensate tends to 
zero. (RHS) Chiral extrapolation of a using the T = 0 extrapolated points from the 
LHS. The nonzero intercept lends support to magnetic catalysis. 
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